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Abstract. We study the absorptive and dispersive properties of Doppler- 
broadened atomic media as a function of detuning. Beginning from the exact 
hneshape calculated for a two-level atom, a series of approximations to the electric 
susceptibility are made. These simplified functions facilitate direct comparison 
between absorption and dispersion, and show that dispersion dominates the atom- 
light interaction far from resonance. The calculated absorption and dispersion are 
compared to experimental data, showing the validity of the approximations. 



PACS numbers: 32.70.Jz, 42.62.Fi, 42.50.Gy 

1. Introduction 

In his paper we investigate the relationship between the absorption and dispersion 
experienced by ofF-resonant radiation interacting with an inhomogeneously broadened 
atomic medium. Of particular interest is the response of the medium when the 
detuning is larger than the inhomogeneous linewidth. In this region scattering of 
photons is reduced, but this does not necessarily mean that the dispersive atom-light 
coupling suffers accordingly. The large dispersion of off-resonant Doppler-broadened 
systems has been exploited in a number of studies. Slow-light experiments [Tl |2] 
and theoretical studies [3] utilise the large group refractive index associated with 
the medium to control the propagation of broadband optical pulses. Slow-light 
interferometers using monochromatic [H [5] and broadband light [6] have also been 
demonstrated. The off-resonant Faraday effect can be used to separate the sidebands 
of Raman light with high fidelity [7], and can be used as a dispersive probe with 
continuous- wave or pulsed light [8] . The non-invasive nature of off-resonant dispersive 
probing could also lead to the possibility of "weak" measurements [9] in the context 
of quantum non-demolition (QND). In many of these experiments there is a trade- 
off between the absorption of the optically-thick medium and the magnitude of the 
effect of interest. As the real and imaginary parts of the medium's susceptibility show 
different spectral dependences it is not obvious at which detuning to perform the 
experiments. The motivation of this study is to take previously developed analytic 
results for the susceptibility [10] and to investigate the domain of validity of two 
approximations which facilitate the analysis of experimental data. 

The structure of this paper is as follows. In section [5] we describe the function 
which governs the absorption and dispersion of a Doppler-broadened medium, and 
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go on to make approximations to this function. We then use these analytic 
approximations to compare the absorptive and dispersive characteristics of an atomic 
resonance. In section [3] we compare the theoretical expressions to experiment, and in 
section!?] we draw our conclusions. 



2. The electric susceptibility 

The electric susceptibility of a medium, describes the medium's absorptive and 
dispersive properties. For the case of an isolated resonance in a Doppler-broadened 
atomic medium the susceptibility as a function of detuning from resonance, A, is jlOj : 

^iA)=cl^^s{A). (1) 

Here c^^^ is the transition strength factor for the transition |i^g,TOFg) ~^ |Fe,TOFo), 
d = (Lg||er||Le) is the reduced dipole matrix element of the \Lg) \Lc) transition, 
M is the atomic number density of state \Fg,mpg), and s(A) is the lineshape factor. 
The total susceptibility of the medium is obtained by summing over all transitions 
which the light is stimulating. The absorption coefRcient is proportional to the 
imaginary part of the susceptibility, x\ a-nd has the form of the well-known Voigt 
profile. Dispersion results from the real part, x^- 

The lineshape factor s(A) is the convolution of /(A), the homogeneous atomic 
lineshape, and g{v), the Gaussian distribution of velocities v. This convolution is given 
by: 

/ + 00 
/(A - kv) X giv)dv, (2) 
-OO 

where k is the wavenumber of the radiation, and 

/(A) ^ (3) 
giv) = J_e-W")^ (4) 

Uy/TT 

Here F is the FWHM (Full-width at half-maximum) of the homogeneous broadening, 
and u is the 1/e of the inhomogenous broadening mechanism (and the RMS atomic 
speed). The lineshape s(A) is related to the Faddeeva (or complex error) function, 
wiiz), of complex argument z, via 

s{A) = ^^«(^z), (5) 
ku 

i /■+°° e"^^ 2 

w(iz) — — dx = e^ erfc(z), (6) 

J-oo iz-x 

Equation ([5]) is the exact analytic lineshape of the Doppler-broadened susceptibility; 
unfortunately this exact result can be difficult to use. Although algorithms exist for the 
Faddeeva and complementary error function erfc(z), they are not easy to manipulate 
analytically, and can be time-consuming to evaluate numerically. Consequently it is 
difficult to relate x to z and the parameters of which they are composed (namely 
the widths F and ku). This in turn makes it difficult to see the relationship between 
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the absorptive and dispersive properties. The preceding reasons motivate our study 
of approximations to the analytic result by looking at the Faddeeva function in two 
regimes, named the Gaussian and Lorentzian approximations for reasons that will 
become apparent. 



2.1. The Gaussian approximation 

We consider the situation where the broadening due to atomic motion is much larger 
than natural broadening, which is the case for typical room temperature alkali-metal 
atoms. For this approximation we therefore look at the limit that T /ku — s- in the 
derivation of the Faddeeva function. Starting from the homogeneous lineshape of 
equation ([S]), 

kii A 

lim f{z)^ -— + i^5{ — ), (8) 
r/fc„^o ' ^ku" ^ ' 

where the imaginary part is given by a Dirac delta function. With this expression 
substituted into equation ([5|), the real and imaginary parts of the susceptibility are, 
respectively 



Y^e-(^/'=")'erfi(A/fcu), (9) 
ku 



= ^e-(^/'=")'. (10) 
ku 

The real part contains the imaginary error function erfi(z) which is similar to the 
Faddeeva function in that it needs to be evaluated numerically. The imaginary term is 
the convolution of a Gaussian and a Dirac delta function, and as expected evaluates to 
the Gaussian function responsible for Doppler-broadening, with the FWHM Doppler 
width Awd = 2-\/ln 2 ku. 

In this approximation both and s^ have a Gaussian detuning dependence, 
whose exponential decrease means that it decays rapidly away from resonance, 
has an additional imaginary error function dependence which increases rapidly with 
detuning; hence dispersion contains the long-range characteristics associated with the 
Faddeeva function. Thus the imaginary part will only be valid close to resonance, 
whereas the real part of the Gaussian approximation is expected to be in good 
agreement with the Faddeeva function for all detunings 



2.2. The Lorentzian approximation 

In the Gaussian approximation we made the assumption that homogeneous broadening 
was negligible compared to inhomogcneous broadening, based on the ratio of their 
frequency widths. We saw that this is not true far from resonance. In this section we 
will find regimes under which homogeneous dominates the susceptibility. We begin 
by noting that the complementary error function can be written in the form of a 
continued fraction [TTl [Hj 

/■°° 2 2e"^^ 
erfc(z) = 2 / e~* dt = ^ . (11) 



For \z\ ^ 1 the continued fraction can be approximated to e ^ /z. This requires 
either of the following conditions to be fulfilled: 
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(?) |A| > AcjD 

(a) r > AwD 

The first condition is that the laser is detuned from resonance further than the Doppler 
width and is essentially a property of the light source; the second is that natural 
broadening dominates over Doppler broadening and is a property of the medium. The 
Doppler width can be reduced by, for example, using cold atoms at sub-milliKelvin 
temperatures. Many experiments, including the ones considered in this paper, are 
conducted with alkali- metal atoms on the D-line at room temperature (or hotter) ; the 
parameters of interest are then F ~ 27r x 5 MHz, and Awd ~ 27r x 0.5 GHz, thus 
T / Auj£, ^ 10^^. Therefore, for the limit \z\ ^ 1 to be valid, it is necessary to be 
detuned far from resonance, |A| 3> Acjd- 

Substituting the approximated erfc(z) into ^ we get the result 

i 1 i 

* ~ fc^I ^ r/2-iA' ^ ' 

Note that this is identical to the case for homogeneous broadening, e.g. an ensemble 
of stationary atoms, or atoms at ultralow temperatures for which Doppler broadening 
is negligible [13] . The real part of the susceptibility gives the dispersion function, and 
the imaginary part the Lorcntzian function; specifically 

"""^^ ' (r/2)2 + A2' ^^^^ 

= (14) 

(r/2)2 + A2- ^ ' 

Furthermore, since |A| Ao^d » F, these relations simplify further to = — 1/A, 
and = r/2A2 respectively. These detuning dependences are discussed further in 
section 12.41 

The physical interpretation of the Lorentzian approximation is that the Gaussian 
lineshape responsible for inhomogeneous broadening decreases exponentially with 
detuning, whereas the homogeneous lineshape decreases much more slowly in the 
wings. Hence the contribution to the overall lineshape far from resonance is 
dominated by the Lorentzian function, and both absorption and dispersion will be 
well approximated. 



2. 3. Validity of the approximations 

Figure [1] shows the lineshape, s, of the Faddeeva function and its Gaussian and 
Lorentzian approximations, for a typical room temperature alkali-metal atomic 
ensemble where the Doppler-broadening is two orders of magnitude larger than natural 
broadening. It can be seen in figure^a) that for |A| < 1.5 x Awd the imaginary part 
of the Faddeeva function is adequately described by the Gaussian approximation, and 
for |A| > 2Acjd the Lorentzian approximation holds. Therefore, close to resonance, 
Doppler-broadening dominates the absorptive interaction; whereas natural broadening 
dominates at large detuning. A similar situation for the real part of the Lorentzian 
aproximation is seen in figure [H^b), i.e. it is valid for |A| > 2Au!d. However, the 
Gaussian approximation is in good agreement with the Faddeeva function over the 
whole spectral range, as predicted in section [2Tl 
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2.4- Comparing absorption and dispersion 

Figure [2] shows the ratio |x^/x^|, calculated using the Faddeeva function and its 
Gaussian and Lorentzian approximations. It shows that the ratio between dispersion 
and absorption continually increases with detuning, with the asymptotic limit 2|A|/r 
from the Lorentzian approximation. Note, however, that dispersion also decreases 
linearly in this limit. Hence, any dispersive effects which require low absorption are 
best performed far from resonance under conditions which increase the atom-light 
interaction e.g. high atomic density [H] or stronger coupling [14| . 

2.5. Hyperfine structure 

We have shown that Doppler-broadening can effectively be ignored for detunings 
|A| > 2/S.Lo-Q. However, this situation is somewhat complicated due to the presence 
of hyperfine structure. For alkali-metal atoms the hyperfine structure is such that 
the ground state splitting, Awhfs, is much larger than the room temperature Doppler 
width. This is not the case for the excited states, which tend to have intervals of 
comparable size to Awd. Hence, in order to calculate x near to the line-centre, each 
individual hyperfine transition needs to be modelled individually, although for some 
purposes excited state splitting can be ignored (for example, on the D2 lines of Rb ^ 
and Cs [2]). Far from line-centre, at detunings larger than the ground state hyperfine 
splitting, it is possible to approximate all hyperfine transitions to a single Lorentzian 
function. By performing this calculation we find that there is a less than 5% error for 
|A| > 3.5 X AtJhfs- 

3. Comparison between theory and experiment 

In order to test experimentally the validity of the approximations to the Faddeeva 
function, the transmission of a probe beam on the Di line of rubidium was recorded. 
The output from an external cavity diode laser at 795 nm was attenuated to be less 
than 1 /xW such that it is in the weak probe limit (see ref. [15]) and passed through 
a 75 mm heated vapour cell, based on the design of [TO]. A solenoid provided the 
heating and magnetic field, when required. The cell contained Rb isotopes according 
to the ratio ^^RbiSSRb of 99:1. 

3.1. Absorption 

The solid curve in figure [31(a) shows the transmission measured at 132°C as a 
function of detuning from the weighted line-centre in units of Doppler width, 
Awd = 27r X 584 MHz. Absorption in the region shown is due to the ^^Rb 
Fg = 2 ^ i^c = 1,2 transitions. Theoretical transmission (dashed curves) was 
calculated using equations ([5]), pO]) and (fT4|) to model each hyperfine transition 
involved. The transmission difference between theory and experiment is shown in 
figure [31Jb). It can be seen that the Faddeeva function agrees with the measured 
data to within the noise level; any discrepancy beyond this is due to the fitting of the 
frequency axis. Both the Gaussian and Lorentzian approximations agree on resonance. 
This is because the transition is optically thick so any variation between the two 
approximations is obscured. At a detuning of about one Doppler width the Lorentzian 
is no longer at zero transmission and only matches with measured data again for 
detunings greater than two Doppler widths. Conversely, the Gaussian approximation 
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matches the experiment up to about 1.5 Doppler widths before differing significantly 
for larger detunings. It was stated in the derivation of the Gaussian and Lorentzian 
approximations that Doppler broadening dominates close to resonance, whilst natural 
broadening dominates far from resonance. It can be seen that at around two Doppler 
widths absorption due to the Gaussian lineshape rapidly decreases, and it is from this 
point that the Lorentzian function becomes the dominant broadening mechanism. 

The situation with the two approximations is similar to that seen in previous 
experiments, where Gaussian fits to data are used when the line-centre is of interest, 
e.g. ref. [T7], and Lorentzian fits for off resonant behaviour, e.g. ref. [2[. However, we 
have derived these lineshapes ab initio and quantified their regime of validity. 

3.2. Dispersion 

The dispersive properties of the medium are easily probed via the Faraday effect f^. 
A rotation of the plane of polarisation arises due to the difference in the refractive 
indices of left and right circularly polarised light. In a heated vapour this rotation can 
be as large as tens of radians over a frequency range of many Doppler widths. Using 
the technique described in ref. [5] we measured the rotation of light polarisation using 
a differencing signal in a balanced polarimeter. Light transmitted through the vapour 
is sent through a polarisation beam splitter and the resulting vertical and horizontal 
polarisations directed to a differencing photodiode, where the intensities Ix and ly 
are subtracted. An important feature of this signal is that the period of oscillation 
is dependent on dispersion, whilst being independent of absorption. From the zero 
crossings one can extract the rotation angle, 6. 

Figure ID^a) shows the differencing signal for an atomic temperature of 112°C 
and applied field of 200 G. The signal is normalised to the maximum intensity, Iq, 
received by one of the photodiodes in the absence of a magnetic field. Also shown is 
the theoretical signal calculated by solving the complete Hamiltonian of the system, 
with a Faddeeva lineshape. There is good agreement between the two curves, the 
main difference being in the amplitude, which is due in part to the differencing 
photodiodes not being perfectly balanced. FigurelDJb) compares the Faddeeva function 
and its approximations to the rotation angle experienced by the linearly polarised 
probe beam. The data points are taken from the zero crossings in figure [U^a). 
Excellent agreement is seen between measured data and both the Faddeeva function 
and Gaussian approximation over the whole spectral range, whilst the Lorentzian 
approximation only holds far from resonance, in agreement with the conclusion of 
sections 12.11 and 12.21 

4. Conclusion 

We have seen in the previous two sections that the Faddeeva function gives a good fit 
to data over the whole frequency range, whereas the Lorentzian approximation is only 
valid far from resonance. The Gaussian approximation, however, shows contrasting 
behaviour in its absorptive and dispersive properties. With regard to absorption it is 
valid close to resonance only, whilst the dispersive properties are accounted for over 
all detuning. The differing nature of the two approximations stems from the fact that 
only in the Lorentzian approximation did we assume that |A| 3> Awd- 

The approximations to the electrical susceptibility we have described facilitate 
(i) the analytic manipulation of x, allowing a comparison between the dispersive and 
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absorptive properties of atomic media; and (ii) ease of computation of properties of 
interest. In particular, we have seen that for a detuning greater than two Doppler- 
broadened Hnewidths the fully analytic Lorentzian approximation is valid. From this 
we have shown that off resonance dispersion increasingly dominates over absorption. 
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Figure 1. Comparision of the Faddeeva function and its approximations. The 
real and imaginary parts of the susceptibiUty lineshape, s, are shown in parts (a) 
and (b), respectively. The horizontal axis is detuning, A, in terms of the Doppler 
width, AajQ. The solid red curve shows s calculated using the Faddeeva function 
(with r/AoJD = 10~^), whilst the Gaussian and Lorentzian approximations are 
shown as dot-dashed green and dashed blue curves, respectively. 
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Figure 2. The relative importance of dispersive and absorptive properties of a 
Doppler-broadened atomic medium. The ratio |x^/x'l ^ single transition is 
shown as a function of detuning, A, in terms of the Doppler width, Ao^d. The solid 
red curve shows x calculated using the Faddeeva function (with F/Acjd = 10~^), 
whilst the Gaussian and Lorentzian approximations to x ^'I'e shown as dot-dashed 
green and dashed blue curves, respectively. 
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Figure 3. Comparison between experiment and theory for the transmission 
of a weak probe beam through a vapour cell (a) Experimental data are shown 
in red, whilst the dashed black curve shows the transmission calculated using 
the Voigt function. The Gaussian and Lorentzian approximations to the Voigt 
function are shown as dot-dashed green and dashed blue curves, respectively, 
(b) The difference in transmission between theoretical and measured data. The 
experimental data were obtained with red-detuned light, but plotted against 
A' = —A. The origin of the detuning axis is from the *''Rb _Fg = 2 ^ Fe = 1 
transition. 
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Figure 4. Comparison between experiment and theory for the polarisation 
rotation of a weak probe beam through a vapour cell, (a) Balanced differencing 
signal. Experimental data are shown in red, whilst the dashed black curve 
shows the theoretical signal calculated using the Faddeeva function, (b) The 
rotation angle of the beam's plane of polarisation: red points are from the 
zero crossings of the measured data, curves are calculated using the Faddeeva 
(black), Gaussian (green) and Lorentzian (blue) functions. The Doppler width is 
Aujo = 2n X 569 MHz. 



